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The 3:1 mean-motion resonance of the planar elliptic restricted three body 
problem (Sun-Jupiter-asteroid) is considered. The double numeric averaging is 
used to obtain the evolutionary equations which describe the long-term behavior 
of the asteroid's argument of pericentre and eccentricity. The existence of the 
adiabatic chaos area in the system's phase space is shown. PACS: 95.10.Ce, 

96.30.Ys 

1. Introduction 

The 3:1 mean-motion resonance in the planar elliptic restricted three-body problem 
(Sun-Jupiter-asteroid) has long attracted considerable attention of specialists [1-16]. 
In order to find secular effects, the equations of motion can be averaged over fast 
variables, namely, over mean longitudes of the asteroid and Jupiter (see, for example, 
[3, 12]). Upon averaging, a nonintegrable system appears which describe the "fast" 
and "slow" components of secular evolution. The "fast" evolution consists in changing 
resonance phase (Delaunay variable) D = A — 3A' , where A and A' are the mean 
longitudes of the asteroid and Jupiter, respectively. The "slow" evolution reveals itself 
in a gradual change of perihelion longitudes of resonance asteroid orbits. 

In [4-6,8,11] a model Hamiltonian system was considered with Hamiltonian TC p 
which was a principal part of the Hamiltonian of a planar elliptic three-body problem 
averaged over A and A' taking the resonance into account. Such a truncation of 
the Hamiltonian is justified at small orbit eccentricities of the asteroid and Jupiter. 
Processes described by this model system are also separated into "fast" and "slow" 
processes. In order to analyze different variants of the " slow" evolution, one can make 
yet another averaging: averaging over fast processes [6,17]. 

In this paper double averaging is used for studying the "slow" evolution without 
restrictions on the orbit eccentricity of an asteroid. The results of this study allow 
one to understand, under which initial conditions the system with Hamiltonian TC p 
correctly describes secular effects at the resonance under discussion. 
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2. Averaging over mean longitudes 



We assume that the semimajor axis of the orbit of Jupiter can be taken as the unit 
length, while the sum of masses of the Sun and Jupiter is the unit mass. The unit time 
is chosen so that the period of revolution of Jupiter around the Sun is equal to 2n . 

We write the equations of motion of the asteroid in the variables 

x, y, L, D, 

where x, y, and L are the elements of the second canonical Poincare system, and they 
are related to osculating elements by the formulas 

x = y 2^/(1 - (j,)a[l - Vl - e 2 ] cosu, (2.1) 



y = — y 2^/ (1 - /j,)a[l - Vl - e 2 ] sincu, 
L = yj(l -fi)a. 

Here, uj, e, and a are the longitude of perihelion, eccentricity, and semimajor axis of 
the asteroid orbit, and /i is the mass of Jupiter. 

The equations of motion have the following canonical form 



dx 


dIC 


dy _ 


dIC 


~dt ~ 


dy ' 


dt ~ 


dx ' 


dL 


dK, 


dD 


dK 


~dt ~ 


3D ' 


~dt 


~ ~dL 



(2.2) 



with the Hamiltonian 

/C = -fi^£-3L-Mfl. (2.3) 
Function R in the expression for /C is defined in the following way: 

P 1 (r, V) 

1\, O , 

|r-r'| r' 3 ' 

where r = r(x, y, L, X(D, A')) and r' = r'(A') are the heliocentric radii vectors of the 
asteroid and Jupiter. 

Formal averaging of the equations of motion consists in substitution of the function 

1 r 2lT 

W(x, y,L,D) = — R(x, y, L, D, A') d\'. (2.4) 

Z7T JO 

for function R in expression (2.3) for Hamiltonian. 

As a result, the equations of motion become autonomous: the mean longitude of 
Jupiter A' = t + X' is eliminated from their right-hand sides. 

A detailed description of the numerical algorithm used to construct averaged equa- 
tions at resonance 3:1 in the general case (the motion of the asteroid is not restricted 
to the orbit of Jupiter) is given in [18,19]. In the planar problem the similar formulas 
can be written in less cumbersome way (see Appendix). 
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3. Fast-slow system describing secular effects in the motion of 

a resonance asteroid 



In motion of the asteroid in resonance 3:1 with Jupiter the value of variable L is 
close to L = l/s/3. In the limiting case \x = the asteroid, moving along the orbit 
with the semimajor axis a res = L\ pa 0.48074 , makes exactly three revolutions around 
the Sun during one revolution of Jupiter. 

Following the general scheme of studying resonance effects in Hamiltonian sys- 
tems [20] we change variable L in system (2.2) averaged over A' for variable $ = 
(Lo — I)/ yZ/i representing the normalized deviation of L from its resonance value and 
introduce a new independent variable r = yfjit . Restricting ourselves to the leading 
terms in the expansion in terms of e — yjji of the right-hand sides of the equations of 
motion in variables x,y,$, and D we get: 

dx dV dy dV 
dr dy ' dr dx ' 

dD d$ dV 

= n<X> = 

dr ' dr dD' 

where 

V(x, y, D) = W(x, y, I , D) , a = ^ = 9^3 . 

If we take as conjugate canonical variables xj^fe and y/y/s, D and $, system 
(3.1) is Hamiltonian with the Hamilton function 

H = — + V(x,y,D). 

In the general case, variables x,y,D, and $ have different rates of variation: 

dD rf$ dx dy 
— — ~ 1 — — ~ s. 

dr 1 dr ' dr ' dr 

We will refer to variables D, $ and x, y as "fast" and "slow" variables, respectively. 

System (3.1) allows one to investigate secular effects in the dynamics of resonance 
asteroids without constraints on the orbit eccentricity value. Evolution of osculating 
elements e and uj with an error 0(e) is described by the relations 

e=^vV + 2/ 2 ) [4Io - (^ 2 + I/ 2 )] , (3.2) 

2tt — arccos — — = , y > 

y x 2 + y 2 

(r- ■ if ^ 0). 

arccos , x , y < 

\Jx 2 + y 2 

where variables x and y change in the region S = {(x,y),x 2 + y 2 < 2I } • 

Taking into account the separation of variables into "fast" and "slow" variables, we 
refer to system (3.1) as a fast-slow (FS) system [20]. 
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4. Principal part TL P of FS-system Hamiltonian 

In [4, 6, 8, 17] the properties of a Hamiltonian system with Hamiltonian 7i p were 
studied. After changing some designations it takes on the form: 

a<& 2 

H p = — + V p (x,y,D). 



Here 



D 



V p (x,y,D) 
V = lim e'~ 2 (V(e'x, e'y, D) - V ) - V 2 . 

e'— >0 

Constants V and V 2 in (4.2) have the following values: 



(4.1) 
(4.2) 



Vn 



bf\a res ) 2 



7T 



K(a res ) « 1.0667, 



V 2 = — 



d a 2 db\ 



(o)' 



da \ 2 da 



0.1422. 



In the above formulas o^(a res ) is the Laplace coefficient, and K(a res ) is the complete 
elliptic integral of the first kind with modulus a res . 

Let us assume that (x,y) G Sw = {(x,y),x 2 + y 2 < C^e' 2 } , where Cw ~ 1 . In 
this case the following estimations are valid: 



dV 




dV p 


12 


dV 


dV 


dD 




3D 


~ e , 




~dD 



-J'l 



(4.3) 



dV 




dV p 


~ e', 


dV 


dV p 















- e' 3 , i 



x,y. 



Taking advantage of relationships (4.3), one can demonstrate that, when independent 
variable r changes in the interval ~ e~ l , the behavior of fast and slow variables in the 
solution to FS-system (3.1) which satisfies the condition (x(t), y(r)) G Sw is described 
by analogous components of the solution with Hamiltonian (4.1) at the same initial 



data with accuracies of 



e V 4 and 



e' 3 , respectively. 



Investigation of the system with Hamiltonian (4.1) made by J. Wisdom can be 
interpreted as an analysis of the properties of solutions with slow variables in Sw in 
the case when Cw ~ 5 . Assuming in what follows Cw = 5 , we will refer to Sw as 
the Wisdom region. 

There is an important difference in the behaviors of V(x,y,D) and V p (x,y,D) 
which should be taken into account in future. The Fourier series of function V(x, y, D) 
includes infinite number of harmonics at any (x,y) G S (point (0,0) can be an 
exception: if eccentricity of Jupiter e' = , then V(0, 0, D) = Vq). The expansion of 
V p (x, y,D) into a Fourier series in terms of D contains a finite number of harmonics 
[4]: 

V p (x, y, D) = B(x, y) + A(x, y) cos(D - R(x, y)). (4.4) 

At points 



Pi = (xi,0), x 1 = e'x u X! w 0.14341 



and 



Pi = (x2,0), x 2 = e'x 2 , x 2 
the value of function A(x, y) is equal to zero. 



2.93415 
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5. Properties of the fast subsystem 



At e = the equations of fast variables coincide with the equations of motion 
of a Hamiltonian system with a single degree of freedom which include x and y as 
parameters: 

dD d<S> dV 

The qualitative behavior of trajectories on the phase portrait of system (5.1) is deter- 
mined by the properties of function V(x, y, D) . 

It turns out that there are important distinctions in fast dynamics at e! > and at 
e' < , where ~ 0.0385 is the minimum eccentricity of the Jupiter orbit admitting 
its intersection with the orbit of the resonance asteroid at /i — . 

In the case e' < e' + (orbits of resonance asteroids do not intersect the Jupiter orbit 
at any e < 1) function V(x,y,D) is limited. If (x, y) ^ 5*i U 5*2 , where Si and S2 
are certain 0(e' 3 ) -neighborhoods of points Pi and P 2 , then the phase portrait of fast 
subsystem (5.1) is topologically equivalent to the phase portrait of a mathematical 
pendulum (Fig. la). In what follows we will designate as £0 the set of (x,y) values 
at which dynamics in the fast subsystem is of "pendulum" type. 

The coefficient of the first harmonic in the Fourier expansion of potential V(x, y, D) 
vanishes at points Pi and P2 . Therefore, when (x, y) are selected in the neighbor- 
hood of these points the qualitative behavior of V(x,y,D) is determined by the next 
harmonics of the Fourier expansion. At (x, y) e Si U S 2 potential V(x,y,D) is char- 
acterized (in the general case) by the existence of local minimum and maximum, which 
leads to a more complicated pattern of motion in the fast subsystem (Fig. lb). The 
boundaries of regions Si and S 2 are curves close to astroid, a typical situation taking 
place when reconstructions of the potential depending on two parameters are studied 



In the case e' > , at appropriate choice of the initial value of variable D , the 
motion of the asteroid along the orbit crossing the Jupiter orbit will be accompa- 
nied by formal collisions resulting in divergence of the integrals determining functions 
V(x,y,D) . Let us designate the set of (x,y) values corresponding to the resonance 
orbits with intersections of the orbit of Jupiter as S3 . Figure lc demonstrates an 
approximate pattern of the phase portrait of the fast subsystem at (x, y) G S3 . The 
partition of S into regions with different types of fast dynamics is shown in Fig. 2. 

Let 

$(r,x,y,h) , D(r,x,y,h) (5.2) 
be a solution to Eq. (5.1) satisfying the condition 

H(x, y, $(t, x, y, h), D(t, x, y, h)) = h, 

in which variable D changes in rotational or oscillating manner: 

D(t + T, x, y, h) = D(t, x, y, h) mod27r 

Here, T(x,y,h) is the period of the solution. We associate this solution with the 
following quantity 

q, rT/cr 

I(x, y,h) = — $ (r, x, y, h) dr, 
Art Jo 
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where the value of a is determined by the type of solution. For rotational solutions 
o" = l and, hence, I(x, y, h) is the action integral. For oscillating solutions a — 2 , 
and the value of I(x,y, h) equals a half of the action integral. 

At e 7^ variables x(t) and y(r) in the right-hand sides of Eqs. (5.1) can be 
interpreted as slowly varying parameters. The quantity I(x, y, h) , coinciding to a 
constant factor with the action integral, will be an adiabatic invariant of system (3.1). 

6. Averaging along solutions to the "fast" subsystem 

Averaging the right-hand sides of the equations for x, y in system (3.1) along 
solutions to fast subsystem (5.2) we get evolution equations describing to the error 
0(e) the changes of variable x, y on the time interval with duration ~ 1/e (or ~l//i 
in terms of original time units): 





Here 



' dV\ 1 fT{x, v ,h) dV 

, oC I T(x, y, h) Jo d( 



C = x,y. 

Construction of phase portraits of system (6.1) is the efficient method of studying 
the evolution of slow variables x,y . At different h the phase portraits can differ in the 
number of equilibrium positions and in the behavior of separatrices. In sections 9,10 
we present numerous examples allowing one to get a general idea about the character 
of reconstructions of the phase portrait when h changes. 

7. Auxiliary functions 

Following [17], let us consider the functions 

H*(x,y) = mmV(x, y,D) , H*(x,y) = m&xV(x,y, D). (7.1) 

The behavior of functions H*(x,y) and H*(x,y) in S is determined by the eccentricity 
e' of the orbit of Jupiter. 

At e' — , the values of the functions at point (x, y) depend exclusively on its 
distance from the origin of coordinates: 

H.(x,y) = H*(r), H*(x,y) = H*(r), 

where r = y/x' 2 + y' 2 . The plots of functions H*(x,y) and H*(x,y) in the case e' = 
are presented in Fig. 3a. At the origin of coordinates x = y = functions and 
H* have the same value: 

H*(0,0) = H*(0,0) = V . 

This value is the global maximum and global minimum for functions and H* , 
respectively. The minimum value of functions is reached at the distance 
0.74420 from the origin of coordinates. Function H* reaches its maximum value at 
the distance r* w 1.05431 . 



6 



The plots of functions if* and H* at e' G (0, e'J can be considered as a result 
of deformation of the plots constructed for the case e' — . The periphery parts of 
the if* plot are elevated in the half-plane x > and go down in the half-plane a; < 
(Fig. 3c). Function if* have a global minimum at point P* = (— r* + O(e'),0) , and 
its saddle point is located at P** = (r* + O(e'),0) . The sides of the if* plot are 
displaced in opposite directions: point P* = (— r* + O(e'),0) becomes a point of the 
global maximum, and the saddle is located at point P** = (r* + O(e'), 0) . At e' — > e* 
the point P* is shifted to the boundary of the region, and the value of if* at this 
point increases infinitely. 

In the Wisdom region Sw 

\0(e' 4 ), (7.2) 

/4\ 



ff* 




= V + e' 2 


if* 


f X 


e'J 


+ v 2 


H* 




= V + e' 2 


H* 


/ X 


e'J 


+ v 2 



where 



+ 0(O, 



= nhn V(x, y, D), H — maxF(x, y, D). 



Taking advantage of the results of studying the properties of functions and H 
presented in [17], one can establish that function has the global maximum at a 
certain point P' 2 e 5*2 , a local minimum at point P{ e Si , and a saddle at the point 

Ph = (xh, 0), = e'x/f + o(e'), a;^ ~ 1.33646. 

Function H* will have the global minimum at point P" G Si , a local maximum 
at point P 2 " G S 2 , and a saddle at the point 

P s = (x s , 0), x s = e'x s + o(e'), ^5 ~ 1.86718. 

We will designate the minimum (maximum) values of functions if* and H* in S 
as h^min and /i^ in (h* max and /i^ ax ), respectively. The value of function if* at point 
P* is ft,** , and the value of H* at point P* is ft,** . 

At e' > e* the behavior of H*(x,y) changes substantially: H*(x,y) = 00 for 
(x,y)ES 3 (Fig. 3c). 

8. Forbidden area and uncertainty curve on phase portraits 

The region 

M(h) = {(x,y)eS,H*(x,y)>h} (8.1) 

is forbidden for phase trajectories of system (6.1). At a given h the slow variables 
cannot assume values from M(h) . 
The curve 

r(h) = {(x,y)eS,H*(x,y) = h} 

is called the uncertainty curve. In the case H*(x,y) = h the trajectory of the fast 
system is a separatrix and, hence, one cannot use averaging. If T(h) is present on the 
phase portrait of system (6.1), it consists of several fragments undergoing a series of 
bifurcations when h is varied. 

When a projection of the phase trajectory of system (6.1) onto the plane x, y 
crosses the curve T(h) , a quasi-random jump of adiabatic invariant I(x,y,h) occurs 
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[7]. When studying the evolution of slow variables on a time interval of order of 1/e , 
this violation of adiabatic invariance is usually neglected, and solutions of the averaged 
system on curve T(h) are glued in accordance with the following rule 

Ibefore I after) 

where Ibefore is the value of I(x, y, h) along the part of the phase trajectory of system 
(6.1) approaching T(h) , and hefore is the value of I(x,y,h) on the trajectory part 
going away from curve T(h) . For most initial conditions, the accuracy of such an 
approximation is 0(e) on the specified time interval. 

The phenomena taking place at multiple intersections of the uncertainty curve will 
be discussed in Sec. 13. 

9. Studying slow evolution based on averaged equations: the 

case e' < 

Figures 4-11 present examples of phase portraits constructed with the help of nu- 
merical integration of Eqs. (6.1) under an assumption that Jupiter moves along the 
orbit with eccentricity e' = 0.02 < e'„ . 

For better visualization the phase portraits present the behavior of quantities x 
and y which are related to variables x, y and osculating orbital elements e, cD as 

x i 

x = —J4L - (x 2 + y 2 ) = ecosu, 
ZLq 

y= -^—JaLq - (x 2 + y 2 ) = -ecoso). 

ZLq 

In the case h G (h* m i n ,h**) system (6.1) has one stable position of equilibrium (Fig. 
4). Phase trajectories are located in the region S/M(h) which shrinks to point P* at 

h > ^*mm • 

Bifurcation at h — h** leads to the appearance on the phase portrait of one more 
position of equilibrium (Fig. 5). The forbidden area M(h) is split into two components 
(central and periphery zones). With increasing h the periphery zone is gradually 
contracted and disappears. 

In the interval Iw = {h* m in ~~ CVe' 2 , h* max + Cw^' 2 ) a series of bifurcations occurs 
in which changing behavior of phase trajectories in Sw , and of geometry of T(h) 
and M(h) approximately corresponds to a scenario described in [17]. Insignificant 
distinctions are caused by the difference (discussed in Sec. 4) in the properties of 
function V(x,y) and of V(x,y) in the neighborhood of points Pi and P 2 . They reveal 
themselves in two very narrow intervals of h values ( ~ e' 4 ) whose lower boundaries 
are minimum values of H* on the boundaries of regions Si and S 2 , respectively. The 
approximate formula 

h = V + e' 2 (h-V 2 ) (9.1) 

allows one to relate bifurcation values of parameter h G Iw with bifurcation values of 
parameter h of the Wisdom system calculated in [17]. 

As an illustration, Fig. 6a presents the phase portrait of system (6.1) at h — 
1.068 G Iw ■ Figure 6b gives its enlarged fragment that demonstrates the behavior 
of phase trajectories in the Wisdom region. One can notice a similarity of Fig. 6b 
with Fig. 5 of [6] and Fig. 11 of [17] where phase portraits of the Wisdom system are 
presented. 
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At h ~ 1.119 a reconnection of separatrices takes place (Figs. 7-9). At h ~ 
1.151 , the pair of equilibrium positions located on the right of the Oy axis (Fig. 10) 
disappears from the phase portrait. At h 1.170, the pair of equilibria to the left 
of this axis merges. After that, a single stable position of equilibrium remains in the 
central part of the S region (Fig. 11). 

Further increase of h results in the appearance of equilibrium positions of the 
periphery of S region (Fig. 12). Changing geometry of the uncertainty curve at 
h = h** and its disappearance from the phase portrait at h — h* max are accompanied 
by merging of the positions of equilibrium at points P** and P* , respectively. 

10. Investigation of slow evolution based on averaged 
equations: the case e' > 

If e' > e' m , one needs to take into account the existence of resonance orbits crossing 
the orbit of Jupiter (the orbits with parameters from region S3 ). Region S3 in variables 
x, y looks like a narrow belt on the periphery of region S . If one chooses (x, y) G S3 in 
system (3.1), there are, in the general case, two "fast" processes over which averaging 
is possible. Therefore, the right-hand sides of Eqs. (6.1) are ambiguously determined 
in S 3 . At numerical integration, in the situation when evolving orbit of a resonance 
asteroid begins to intersect the orbit of Jupiter, one can choose an appropriate solution 
to fast subsystem (5.1) as a closest to the solution along which averaging was made on 
the preceding step. 

Figures 13-16 present the phase portraits of system (6.1) constructed for the case 
e' = 0.048 . To choose such value of eccentricity of Jupiter is traditional for numerical 
investigations of the dynamics of asteroids in the context of the restricted elliptical 
three-body problem [4]. 

The magnified fragments of the phase portrait at h — 1.07146 (Figs. 13b and 13c) 
resemble the phase portraits presented on Figs. 9 and 10 of paper [17]. These figures 
show that system (6.1) has the families of closely located equilibrium positions. 

Figures 14-16 demonstrate the reconnection of separatrices at h 1.0955. They 
resemble Figs. 7-9 illustrating a similar bifurcation in the case when Jupiter moves 
along the orbit with eccentricity e' = 0.02 . Nevertheless, there are important qualita- 
tive distinctions: uncertainty curve T(h) intersects seperatrices in Figs. 14-16, while 
there is nothing of this kind in Figs. 7-9. The location of the uncertainty curve on the 
phase portrait determines the size of the region of adiabatic chaos in the phase space 
of a nonaveraged system (Sec. 13). 



11. Bifurcation diagrams 

Numerical integration of Eqs. (6.1) shows that in most cases qualitative changes in 
the behavior of trajectories on phase portraits are caused by some changes in properties 
of the equilibrium solutions lying on the Ox axis. Therefore, when studying the 
bifurcations of phase portraits, of fundamental importance is the analysis of diagrams 
that demonstrate the types and position of equilibrium solutions on the Ox axis at 
different values of h . 

As an example, Figs. 17a and 18a present the bifurcation diagrams constructed 
under an assumption that the eccentricity of Jupiter e' was equal to 0.02 and 0.048 , 
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respectively. The left parts of these diagrams contain the families of stable equilibrium 
solutions generated for h = h* min at point P* and the families of unstable equilibrium 
solutions generated for h = h** at point P** . With increasing h these equilibrium 
solutions are merged with equilibrium solution of two other families originating in 
the central part of S . The right-hand parts of the diagrams indicate to merging 
of equilibrium solutions that takes place at critical points of function H*(x,y) . For 
clearness Figs. 17b and 17c present enlarged fragments of the bifurcation diagram in 
the case e' = 0.02 . In the case e' = 0.48 there is only one similar bifurcation in the 
right-hand part of the diagram: at h = h** . 

The segments of bifurcation diagrams describing reconstructions of phase portraits 
in the Wisdom region are practically merged in Figs. 17a and 18a when h varies in the 
interval Iw ■ A blowup image of corresponding fragment of the diagram for e' = 0.048 
can be seen in Fig. 18b. If one compares Fig. 18b with a bifurcation diagram of the 
Wisdom system presented in [17], some distinctions (along with doubtless similarity) 
should be noted. In particular, Fig. 18b indicates that the stability of equilibrium 
solutions in the lowest family changes at h ~ 1.0716 , i.e., when the values of function 
H* at point P s are exceeded insignificantly (H*(x s ,0) ~ 1.0715). In the Wisdom 
system, a similar bifurcation occurs at the parameter value substantially different from 
the value of auxiliary function at the saddle point. It is quite likely that the differences 
are due to the presence of singular points of function V at e! > . In their neighbor- 
hood the value of this function increases infinitely. The singularities can substantially 
reduce the size of the region where the function's behavior can be studied based on 
consideration of its principal part. 

Further detailed elaboration of the bifurcation diagram description is related to 
the analysis of reconstructions of phase portraits in the Si region at h G Ii — 



{H*{xi,0),h* min + O(e' 4 )) and in S 2 at h G I 2 = (K max , H*(x 2 , 0) + 0(e' 4 )) . It is 



more complicated to seek equilibrium solutions in these cases: as in the S 3 region, the 
right-hand sides of averaged equations can be determined ambiguously. We think that 
such a situation in the fast-slow system of form (3.1) deserves special consideration. 
Therefore, we do not present here the results of performed numerical investigation. 
Since intervals ii and I 2 are very small (~ e' 4 ), the absence of information on bi- 
furcations in Si and S 2 will be insignificant for a specialist who is going to use the 
results of this work when studying the dynamics of resonance asteroids. Even in the 
enlarged fragment of the bifurcation diagram in Fig. 18b one cannot see any families 
of equilibrium solutions in Si and S 2 at h G ii and h G I 2 , respectively. 

12. Periodic solutions and invariant tori of the three-body 
problem at the 3:1 mean motion resonance 

In the general case, the stable (unstable) periodic solutions to FS-system (3.1) and 
stable (unstable) invariant tori in the extended phase space x,y, L, D, \' of original 
system (2.2) correspond to the stable (unstable) stationary solutions to evolution equa- 
tions (6.1). 

Taking advantage of the fact that points 



P* (•£*) 0), P^ 



(x**,0), P s = (x s ,0) 



and 



P* = (x*,0), P 



(XH,0) 
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are critical points of functions if*(x, y) and H*(x, y) , respectively, 
existence of the following stationary solutions to F'S-svstem f3 1^ 

t^vVlO U^IJ.^^ \JL 1U11UVV 111k O U C1j U IWI-IOj-I V OUlUllUllO UW 1 ■ .J OV OllSlll I U . X J 


one can prove the 


x = x*, y = 0, D = n, $ = (h = /i* min ), 


(12.1) 


x = x**, y = 0, D = tt, $ = (/i = /i**), 


(12.2) 


x = xh, y = 0, D = 0, $ = (/i = /in), 


(12.3) 


z^z*, y = 0, L> = 0, $ = 0(/i = /i; a J, 


(12.4) 


x = x**, y = 0, D = 0, $ = (/i = h**), 


(12.5) 


x = x s , y = 0, D = tt, $ = (h = h s ), 


(12.6) 


in which the values of variables x, y form a stationary solution to evolution equation 



(6.1). 

Periodic solutions of the restricted planar three-body problem correspond to sta- 
tionary solutions (12.1)— (12.6). Solutions (12.3) and (12.6) are similar to periodic 
solutions studied by Hill [1] and Sinclair [2]: in these solutions an asteroid moves in 
an nearly circular orbit with eccentricity e ~ e' . In the remaining cases the asteroid 
motion proceeds in orbits with eccentricity e ~ 1 (Fig. 19). The properties of such 
orbits were examined in [10, 13, 14]. 

13. Adiabatic chaos 

In the neighborhood of uncertainty curve T(h) the projection of a phase point of 
system (3.1) onto the plane x, y jumps from one trajectory of averaged system (6.1) 
to another in a quasi-random way: \I after — hefore\ ~ £ • As a result of a series of such 
jumps, the phase trajectories of system (3.1) with close initial data 

|A(0) - D 2 (0)\ w e, |$x(0) - $ 2 (0)| w e, \ Xl {0) - x 2 (0)| w e, \ yi (0) - y 2 (0)| w e, 

can go away to a distance of ~ 1 . Their projections onto the plane x, y will fill the 
region S(/i) which is a set of all trajectories of evolution equations (6.1) intersecting 
r(/i) (Fig. 20). In the phase space of FS-system (3.1) diverging trajectories will be 
located in the region 

£*(/i) = {x, y,D,<Z> : H(x, y, D, $) = /i, (x, y) G S(/i)} . 

We call S*(/i) the region of adiabatic chaos: the complex behavior of trajectories 
in S*(/i) is associated with violation of adiabatic invariance in the neighborhood of 
T(h). 

The properties of adiabatic chaos in Hamiltonian systems were studied in [17, 21]. 
The existence in this region of numerous (~ 1/e) stable periodic solutions was proved 
in [21]. In [17], using numerical methods, such solutions were sought for an autonomous 
FS-system with two degrees of freedom. 

The diverging trajectories of original (unaveraged) system (2.2) correspond to the 
trajectories of system (3.1) diverging in S(/i) . Thus, the regions of chaotic dynamics 
originating due to violations of adiabatic invariance will also exist in the phase space of 
system (2.2). The stable periodic solutions to FS-system (3.1) turn into stable invariant 
tori in the extended phase space x, y, L, D, A' . 
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Conclusions 



Studies of 3:1 mean- motions resonance are of great importance for understanding 
the evolution of orbits of many celestial bodies. Asteroids of the Hestia family move in 
resonance 3:1 with Jupiter. In [22], the possibility of such a resonance was considered 
for Uranus' moons Miranda and Umbriel. A hypothesis of planet motion in resonance 
3 : 1 in the system 55 Cancri was discussed in [23] . 

The approach used in this paper allows us to get a sufficiently detailed description 
of secular effects in motion of resonance objects in the context of a planar restricted 
elliptical three-body problem. 
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Appendix: calculation of derivatives of function W 



1. Preliminary remarks.For calculation of W it more convenient to use instead of 
(2.3) the relation 

W =JT * R(x,y,L,D,\'(\,D))d\. (Al) 

07T JO 

Following [18, 19], let us make in (A.l) a change for integration over eccentric longitude 
g equal to a sum of the eccentric anomaly and the longitude of pericenter of an asteroid 
orbit 

W 1 f 67T R 9X d 
Git Jo dg 

When integrating over g there is no necessity to solve the Kepler's equation for seeking 
the true longitude characterizing the position of an asteroid. The formula connecting 
the mean longitude A with g has the following form: 

A = g — esin(g — to). (A2) 

Calculation of the derivatives of W with respect to x, y, L, and D is reduced to 
differentiation of the integral over parameter 

(dRdX | R d 2 X \ d 
\ dx dg dgdx J 

(dRdX +R d 2 X\ ^ 
\dy dg dgdy ) 



dW 


= -f 


dx 


67T JO 


dW 


= -f 


dy 


6"7T JO 


dW 


= -/' 


dL 


6n Jo 



6 TV 



'dRdX d 2 X \ 
dLdg- + R dgjL) 



dW J_ r** (dRdX\ 

~6Wo \dDd~q) 9 - 



dD 6n Jo \ dD dg t 

2. Derivatives of R with respect to x, y,L,D . In what follows we use the heliocentric 
coordinate system 0^r]( (all bodies move in the plane 0£r) , and the orientation of 
axes is kept invariable). If projections of r onto axes 0£ and Or] are designated as 
£ and rj and projections of r' as £' and rf , then 

dR_dRd^_ dRdji dRd^_ dRdrf_ 
dx <9£ dx drj dx dC,' dx drj' dx ' 

dR _ dRd£ dRdji dRd£_ dRdrf_ 
dy dt; dy drj dy d£' dy drj' dy ' 
dR_dRdi_ dRdr]_ dRd^_ dRdrf_ 
~dL ~ ~dJdL + 1h)~dL + ~d^'~dL + d^'dL' 

dR _ dR d? dR drf 
dD ~ d^dD + d^fdD' 

3. Derivatives of £, rj with respect to x,y,L . One can write the derivatives of £ 
with respect to x, y, and L that appear in (A. 3) in the following form: 

<9£ <9£ de <9£ du 

dx de dx du dx 
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<9£ <9£ de d£ doo 
dy de dy duo dy 
dt; d£ da d£ de 
dL da dL de dx 

Expressions for derivatives of rj with respect to x, y, and L have a similar structure. 

4- Derivatives of £ and rj with respect to Keplerian elements. The coordinates of 
an asteroid £, rj , Keplerian elements of its orbit a,e,cu , and eccentric longitude g are 
connected by the relationships 



£ = a cos g — ecosuo + (1 — Vl — e 2 ) sin(g — a;) sincu 



rj = a smj-esmw- 
After simple transformations we get: 
<9£ £ <9£ 



<9a a' de 

di = 

duo 

drj rj drj 
da a ' <9e 

drj 

duo 



= —a cos uo 



1 — e 2 ) sin((7 — c<j) cos a; 



: sin(g — uo) sin a; , 



VT 3 ^ 



e sina; + (1 — Vl — e 2 ) sin(g — 2uo) 



= —a sincu + 



: sin(g — uo) cos uo , 



e cos uo — (1 — Vl — e 2 ) cos(g — 2cj) 



5. Derivatives of elements of the asteroid orbit with respect to x,y,L . The formulas 
relating x,y,L , and Keplerian elements uo, a, e, have the following form: 



CO 



a = 



— arctg (|Q , x > 
7r — arctg , x < 

l 2 y^ 2 + y 2 ) 



(AA) 



e = 



where 



1 - /j, ' ~ 2L 
s = 4L - x 2 - y 2 . 



Using formulas (A. 4) we find: 
(9a; y 



duo 



x 



da 



2L 



dx x 2 + y 2 ' dy x 2 + y 2 ' dL 1 — /i ' 



(A5) 



<9e 



x 



<9x 2L I y x 2 + y 2 



I x 2 + y 2 



'x 1 + y 1 



de y 

dy 2L \ \ x 1 ■ y 2 v 

<9e \/x 2 + y 2 I 1 V^ N 



<9L 



V^ 2L 
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6. Derivatives of £',r/ with respect to x,y,L,D. These derivatives are calculated 
according to the formulas 

d£' d£' dg' drj drf dg' 
dx dg' dx ' dx dg' dx 

where g' is the eccentric longitude of Jupiter. 

Differentiating the expressions for coordinates of Jupiter 

£' = cos g' — e', rf = \J\ — e' 2 sin g', 

we get: 

— = - smg , — = VI - e 2 cosg . 
We find the value of g' from the Kepler's equation 

g — e' sing' = A' . 

Here, 

A' = ±(A - D) = 1 [g - esm( 5 — ui) — D\. 
Derivatives of g' with respect to e, uj, D have the form: 

dg' sin(g — u3) dg' ecos(g-uj) 



de 3(1 — e'cos g') ' du 3(1 — e'cosg') ' 
dg' 1 



&D " 3(1 - e'cos g') ' 
The values of derivatives of g' with respect to we get x,y,L using formulas (A. 5): 

dg' dg' de dg' du dg' dg' de dg' dtu dg' dg' de 
dx de dx du dx ' dy de dy du dy ' dL de dL 

7. Derivatives of the mean longitude of an asteroid. Differentiating relation (A. 2) 
with respect to eccentric longitude, we find: 

— — = 1 — e cos(o — cu) — 1 — ^—{x cos q — y sin q). 
dg y ' 2L y ' 



Accordingly: 



d 2 X x y/s 



x cos g — y sin g) — — — cos g 



dgdx 2L^s~ K y tf 2L 

d 2 \ y Js 

[x cos g — y sin g) + —— cos g , 



dgdy 2L^ y * y ^ 2L 

d 2 \ 1 / 1 Vi\ 

-- -p - — ](x cos # - y sm g) 



dgdL 2L 
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Fig. 1. Phase portraits of a fast subsystem: 
(a) - (x,y) e S , (b) - (x,y) eS 1 US 2: (c) - (x,y) E S 3 . 
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Fig. 2. Partitioning of set S into regions of the values of slow variables x, y 
corresponding to the dynamics of qualitatively different types in the " fast" subsystem: 
(a) - general view, (b) - enlarged fragment containing region Si . Region S2 has similar 
size (~ e' 3 ) and shape. 




Fig. 3. 3D graphs of the functions H*{x,y) and H*(x,y): 

(a) e' = 




Fig. 3. 3D graphs of the functions H*{x,y) and H*(x,y): 

(b) e' = 0.02 




Fig. 3. 3D graphs of the functions H*(x,y) and H*(x,y) 

(c) e' = 0.048 




Fig. 4. Phase portrait of system (6.1) at h e (/i* m m, ^**), e ' £ (0, e'J . 
Here and below, forbidden area M(h) is shown by black color. 




Fig. 5. Phase portrait of system (6.1) at h e (/i**, hn), e' G (0, e'J . 





Fig. 6. An example of the phase portrait of system (6.1) at h e Iw, e' G (0, e'J 
(a) and an enlarged fragment of this phase portrait demonstrating the behavior of 
trajectories in the Wisdom region S w (b). Thick line here and below represents 
indeterminacy curve T(h) . The region of the values of slow variables at which the 
motion in fast subsystem (5.1) has rotational character at the level H = h is marked 
in gray. 




Fig. 7. Phase portrait of system (6.1) after a series of bifurcations occurring in region 

Sw when h increases. 




Fig. 8. Phase portrait of system (6.1) at bifurcation value of h: there are asymp- 
totic solutions with the different limits at r — > — oo and at r — > +00 . Such solutions 
disappear at arbitrarily small changes of h . 




9. Phase portrait of system (6.1) after reconnection of separatrices. 




Fig. 10. Phase portrait of system (6.1) after merging of equilibrium positions in the 

right-hand half-plane, e' G (0, e'J . 




Fig. 11. Phase portrait of system (6.1) after merging of equilibrium positions in the 

left-hand half-plane, e' G (0, e'J . 




Fig. 13. Phase portrait of system (6.1) at h E I w (e' > e'J (a) and its enlarged 
fragments (b and c). Dashed line shows the position of the boundary between regions 
So and S 3 . 




Fig. 14. Phase portrait of system (6.1) before reconnection of separatrices (e' > e'J 




Fig. 15. Phase portrait of system (6.1) for bifurcation value of h (e' > e'J 




16. Phase portrait of system (6.1) after reconnection of separatrices (e' > 




Fig. 17. Bifurcation diagram in the case e' = 0.02 < e' + (a) and its fragments (b 
and c). For different values of h the diagram presents equilibrium solutions to system 
(6.1) located on the Ox axis. Thin solid and dashed lines correspond to the families of 
stable and unstable solutions, respectively. Thick line shown location of the points at 
which indeterminacy curve T(h) intersects axis Ox. Vertical segments of the region 
marked in black are the segments of axis Ox lying in the forbidden area M(h) at a 
corresponding value of h . 
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Fig. 18. Bifurcation diagram in the case e' = 0.048 > (a) and its enlarged 
fragment corresponding to variation of h in the interval Iw (b). 
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Fig. 19. The value of eccentricity e in periodic solutions to the three-body prob- 
lem corresponding to equilibrium solutions to system (3.1). Solid and dashed lines 
correspond to the families of stable and unstable solutions, respectively. 



